Spin and charge transport induced by gauge fields in a ferromagnet 
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We present a microscopic theory of spin-dependent motive force ("spin motive force") induced 
by magnetization dynamics in a conducting ferromagnet, by taking account of spin relaxation of 
conduction electrons. The theory is developed by calculating spin and charge transport driven by 
two kinds of gauge fields; one is the ordinary electromagnetic field ^Ji™' ''■'^'^ the other is the effective 
gauge field A^^ induced by dynamical magnetic texture. The latter acts in the spin channel and gives 
rise to a spin motive force. It is found that the current induced as a linear response to j4^ is not 
gauge-invariant in the presence of spin-flip processes. This fact is intimately related to the non- 
conservation of spin via Onsager reciprocity, so is robust, but indicates a theoretical inconsistency. 
This problem is resolved by considering the time dependence of spin-relaxation source terms in 
the "rotated frame", as in the previous study on Gilbert damping [,J. Phys. Soc. Jpn. 76, 063710 
(2007)]. This effect restores the gauge invariance while keeping spin non-conservation. It also gives 
a dissipative spin motive force expected as a reciprocal to the dissipative spin torque ("/?-term"). 
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I. INTRODUCTION 



Manipulation of magnetization by electric currentsi"— 
has been studied intensively for a decade because of 
promising spintronic applications)^ Among them, it was 
demonstrated theoreticalljs^ and experimentally^ that 
an electric current in a conducting ferromagnet can 
drive magnetic textures such as domain walls and vor- 
tices. This is understood as due to spin torques that a 
current exerts on magnetization through a microscopic 
exchange interaction. They include the spin-transfer 
torquefi"— which is based on the conservation of total an- 
gular momentum, and its dissipative correction called j3- 



term 



11-17 



which arises in the presence of spin-relaxation 
processes in the electron system. 

In 1986, Berger predicted a reciprocal effect that a 
moving domain wall accompanied by a periodic rota- 
tion of magnetization generates an electromotive force, in 
analogy with the Josephson effect of superconductivity J^ 
This effect is now understood as a motive force acting in 
spin channel, hence called spin motive force^^Ti^ which 
drives majority-spin and minority-spin electrons in mu- 
tually opposite directions. It is also understood to arise 
from a time-dependent magnetic texture in general. Re- 
cently, it was experimentally detected by Yang el al^ for 
a vortex wall in a ferromagnetic nanowire. Similar phe- 
nomena have also been studied in systems with interface 
or nanoparticlesj^"— 

A theoretical framework for studying spin motive force 
in ferromagnets was presented by Volovik,'l2, or earlier by 
Korenmann et al^ To treat electrons in a spin (or mag- 
netization) texture, they introduced a local spin frame 
whose quantization axis coincides with the local spin 
direction,— n; then there arises naturally an effective 
U(l) gauge field, A^, acting in electron's spin channel. 



which gives rise to an effective 'electric' fieldi^ii 



19.24 



E^. 



-Jd.A^o 



doAI)^—n-id,nxh), (1) 
2e 



or a spin motive force, Fg = —eEg (— e : electron charge). 
Recently, it was pointed out that it acquires a dissipative 
correction22i^ 



E' 



= P—h-diTi, 
2e 



(2) 



in the presence of spin relaxation of conduction electrons. 
The total field is then given by E^ = E° + Ef"". These 
two terms are reciprocals to the spin-transfer torque and 
the spin torque /3-term, respectively)^""— and the di- 
mensionless parameter /3 is the same as that of spin 
torque fiir— 

A spin motive field Eg induces an electric current 



j = a^Es + (Ji{-Es) = asEs, 



(3) 



where a-^ (crj^) is a conductivity of majority- (minority-) 
spin electrons, and cTg = cTf ^ o"^ is the 'spin conductiv- 
ity'. In most theoretical studies, this relation is used to 
identify a spin motive force<^""— In the presence of spin- 
orbit coupling, it induces in addition a charge Hall cur- 
rent, asan x E^, where ctsh is a spin Hall conductivity^^ 
and as a reciprocal to this, a spin Hall current induced by 
external electric field will exert a spin-transfer torquei^ 
Enhancement of magnetization damping due to induced 
spin current was also discussed i^^i'^^ 

The purpose of this paper is to develop a microscopic 
theory of spin motive force basing on the gauge field 
mentioned above. For this, we found it instructive to 
treat spin and charge channels in parallel. We thus study 
spin and charge transport induced by two kinds of gauge 
fields, one acting in charge channel (ordinary electromag- 
netic field) and the other acting in spin channel (spin 



motive field). Particular attention is paid to the effects 
of spin relaxation of conduction electrons. 

In the first part of this paper, we study spin and charge 
transport in a uniformly magnetized state induced by an 
ordinary electromagnetic field. Our calculation is equiva- 
lent to the well-studied two-current model,—"— but some 
interesting crossover is pointed out in diffusion modes. 

In the second part, we study a spin motive force by 
calculating electric and spin currents induced by mag- 
netization dynamics. We encounter a difficulty that the 
current induced as a linear response to the effective gauge 
field Af^ contains gauge non-invariant terms in the pres- 
ence of spin-ffip processes. This difficulty is resolved by 
noting that there is another contribution from the source 
term of spin relaxation, as realized in the study of Gilbert 
damping Ji We also found that such additional contribu- 
tion reproduces the dissipative spin motive force. 

Such additional contributions may look tricky, but 
their necessity can be understood on general grounds. 
In the present gauge-field formalism, in which spin and 
charge channels are treated equally, spin conservation 
and gauge invariance (in the spin channel) are equiva- 
lent at the linear-response level because of Onsager reci- 
procity. However, the former is violated by spin-fiip pro- 
cesses whereas the latter should always hold in order for 
the theory to be consistent. These contradictory aspects 
can only be reconciled by some additional contributions. 

The paper is organized as follows. After describing a 
model in Sec. II, we examine in Sec. Ill the density and 
current response to the ordinary electromagnetic field, 
A™. Here the magnetization is assumed to be static 
and uniform. In Sec. IV, we consider the case that the 
magnetization varies in space and time. By introducing 
another gauge field, Afj, which expresses the effects of 
magnetic texture and dynamics, we examine the density 
and current within the linear response to Af^, with an 
unpleasant, gauge-dependent result. This problem is re- 
solved in Sec. V, where a dissipative correction to spin 
motive force is also obtained. Results and discussion are 
given in Sec. VI, and summary is given in Sec. VII. Cal- 
culational details are given in Appendices. 



where ^^{x) = {cl{x) , c[{x)) is the electron creation op- 
erator at X = {t, r), ef is the Fermi energy, M is the s-d 
exchange coupling constant, n is a unit vector represent- 
ing the direction of d-spiu)^ and <t is a vector of Pauli 
spin matrices. The impurity potential is modeled by 



where U[ and Ri are the strength and position of normal 
impurities, which introduce momentum relaxation pro- 
cesses, and Us and R', are those of quenched magnetic 
impurities with spin Sj, which introduce spin- relaxation 
processes . ^^'^^ We take a quenched average for the impu- 
rity spin direction as 5f = and^ 



S^'S^ 



S,^,S''0 X <! ^ 



si {a,P = x,y) 
S'i {a, l3 ^ z) 



(7) 



as well as for the impurity positions, R'^ and i?'-. When 
the magnetization is uniform and static, n = z, the 
impurity-averaged Green's function is given by 



Gkaiz) 



1 



z- ek + eFa + ij„sgn{lmz) ' 



(8) 



where fe is a wavevector, e^ = h^k^/2m, and epo- = 
ep + ctM. The subscript a =t, i represents the major- 
ity and minority spins, respectively, and corresponds to 
a = +1, —1 in the formula (and to a =i,t or —1, +1). 
Treating Vlmp as perturbation, the damping rate ^^ is 
evaluated in the first Born approximation as 



Jcr = = 7r(rii^o. + T2l^a), 

iTrr 



(9) 



where I'o- = mk-p^/2Tr h is the density of states at epa- 
with fcpcr = \/2rneF^/h and 



Fi = TijU^ -I- UsulSl, 
f 2 = 2n,ul^^, 



(10) 
(11) 



II. MODEL 

We consider a ferromagnetic conductor consisting of 
conducting s-electrons and localized d-spins. We assume 
that the s-electrons are degenerate free electrons subject 
to impurity scattering, and localized d-spins are classi- 
cal, which are mutually coupled via the s-d exchange 
interaction. The Lagrangian for s-electrons is given by 
L = Lc\ — Hsd- 



ifi— + -— V^ + ep - Mmp 
ot 2m 



dr c^ 
Hsd = -M dr n- {c^crc)^, 



c, (4) 
(5) 



with Ui and Ug being the concentration of normal and 
magnetic impurities, respectively. The first and second 
terms in Eq. ^ come from spin-conserving and spin-flip 
scattering processes, respectively. 

In this paper, we assume j^ <^ s-pa and focus on diffu- 
sive transport induced by slowly- varying external pertur- 
bations (electromagnetic fields or time-dependent mag- 
netic texture) . Let q and uj be wavenumber and frequency 
of the perturbation, and define 



Xcr = {D^q'^ - iuj)Ta, 



(12) 



with a diffusion constant D^. Then our assumption 
throughout the paper is expressed as 7^ ^ e^a and 

\XA « 1. 



III. SPIN AND CHARGE TRANSPORT IN 
UNIFORMLY MAGNETIZED STATE 

A. Linear response to electromagnetic field 

Let us examine the density and current response in 
the charge channel, j^ = {p,j), and spin channel, 
js,fi = {ps,js), to the external electromagnetic field, 
AJ;™ = (-<^<=™, A™)4iii^ Here (j)""^ and ^"" are scalar 
and vector potentials, respectively, and the time and 
space components of the four currents are given by 



3=3^" 



ectc(=,r), 
e 



?7i Zmt 



- ^°)^ 



Ps = -ec^<y'ci^j^), 



3s-3i'^ + -PsA--. j(") 



2mi 



(13) 

c^Vc, (14) 

(15) 
c^a'^Vc, (16) 



with c'^^c = c^'^c— {\'c')c. We have defined ps and jg to 
have the same dimensions as p and j, respectively. The 
coupling to the external fields is given by 



-'^em 



dr (p<?i""-j("'-A'= 



(17) 



The currents, j^ and jg.^, are evaluated in the linear 
response to A°™ as 

{j,{q))^ = e2x-(g,c. + zO)A^-(c.), (18) 
(JsAq))^ = e^K^;AQ,^ + ^0)AZA^), (19) 

where A™ (oj) is a Fourier component of A™(x) . The 
response functions X^™ and K^^ are obtained from 



e'KfZiq,iu^x) 



l/T 



dre-^-(T.j(")(g,r)j(")(-g)) 



+ —{p)Sf,^{l-S^o), 



(20) 



e'Klliq^iLo^) 



l/T 



dre-^-(T.jW(g,r)j(°)(-g)) 



+ —(/Os) (5^1/(1 -(^ivo), 



(21) 



by the analytic continuation, iuJx —^ fku + iO, where lox = 
2tt\T (A : integer) is a bosonic Matsubara frequency. In 
this paper, we focus on absolute zero, T = 0. The average 
(• • • ) is taken in the equilibrium state determined by L. 
The Fourier components of the currents are given by 



J^C?) = -e^v^4_^^Cfc_^,<,, 



h.a 



h.a 



(22) 
(23) 



(a)ir 



^W 



(b)A^ 





V, 




Fi f bv^ + r2 f o'Wiy 



Fi ^ 9^1^ + r2 1^ 



FIG. 1: (a) Diagrammatic expression of /CJ;^ . The thick (thin) 
solid line represents an electron line carrying Matsubara fre- 
quency ie„ + ii-o\ {ie„). The shaded part represents the vertex 
function, AJ^. (b) Dyson equation for AJ^. The dotted lines 
represent impurity scattering, either witli (r2) or without (Fi) 
spin-flip scattering. 



with 



1 (m = 0) 

hki/m (/i = i = 1, 2, 3) 



(24) 



and k± = k± q/2. 

The response functions are evaluated with the ladder- 
type vertex corrections^ [Fig.[lja)]. Deferring the details 
to Appendix A, we give the results in the next subsection. 
The results are concisely expressed with the quantities 



Ya = Daq - iu!, 

Z = Y^Y:, + 27:f2{Yiy), 



(25) 
(26) 



and a notation, (••■), meaning to sum over a =t,i; for 
example, (v) = v-^ + h'l, {ah') = i/^ — ly^, {Di') = D^u^ + 
D^Vi, and {oDv) = D^v-^ - D^i/;. By defining {Y)„ = 
Yg., we may also use {DvY) ~ D^v^Y]^ + D^v^Y^ and 
{oDuY) = D^v^Y^ - D^v^Y^. 



B. Result 

1. Charge channel 

The response functions Kf^1{q,uj + iO) [Eq. (PD|) ] for 
the electric density/current are obtained as 



k: 



q^K, 



iO 



k: 



K?, 



where 



lOJ 



K = 



0, - qiUjK, 
(Du) [ 5, 



mi 

n'2 



iujK 



qiqj 



{DiyY)+2TTr2{iy){Diy) 
Y^Y^ + 2nt2{Yu) 



(27) 
(28) 

(29) 



(30) 



The following properties are seen. 

(i) Gauge invariance^ and charge conservation are sat- 
isfied, 



with 



/v-g, = 0, q^K'^l^O, 



(31) 



where g^ = {—to, q) is a four wavevector4^ 

(ii) For r2 = (without spin-flip scattering), we have 






DnVr, 



D^q — i(^ 



(32) 



This means that up- and down-spin electrons diffuse in- 
dependently, and there are two independent diffusion 
modes. 

(iii) For r2 7^ 0, and in the long-wavelength and low- 
frequency limit, T^^ = 2'itT2{v) /h ^ |Yct|, we have 



K = 



{u){Du) 



^c/e^ 



where 



{Dv)q'^ — iuj{iy) D^ffq^ — iuj 



Detf = 



) l^t + n 

is the effective diffusion constant, and 



ac^e^Diy) ^e'Y.Do 



(33) 



(34) 



(35) 



is the electrical conductivity. There is only one diffusion 
mode owing to the spin mixing r2- In the opposite limit, 
Tg7^ ^ I^ctI, we have the behavior ([32]). 

Finally, the charge density p = {jo{q))uj and the cur- 
rent density j = {ji{q))u> are given by 



P 
3 



<7c.E 



2{D'^vY)+2TTV2{Dvf 
Y^Y^ + 2T,f2{Yv) 



(36) 
V(div£;), (37) 



where E{q, u) is a Fourier component of the electric field: 
E{q, uj) = ~iq<l)"''{q, uj) + iujA''"'{q, w) with div£; = iq ■ 
E and V(div£;) = iq{iq ■ E). 



K"" = 



{(jDi^Y) +2'kV2{v){cjDv) 



The difference 

AX" = 27rf2 



Y^Y^ + 2^V2{Yv) 
{aDvY) +2TiV2{av){Dv) 



Y^Y^ + 2^V2{Yu) 
[crv){Dv) - {v){(jDv) 



(42) 
(43) 



Y^Y^ + 2^T2{Yv) 
= 27rf2(crci^- - (7s'^+)/Ze'^ 
= 2TTr2iy+<Jc{P.-Pj)/Ze^ 

arises if f2 7^ (and P^ ^ Pj)- In Eq. (glD, 



e^ ^ (iDcriya 



(44) 



(45) 



is the 'spin conductivity', and P^ = v^Jvj^ and Pj = 
tTs/tTc represent spin asymmetry in the density of states 
and in current density, respectively, which are different 
in general. The following properties are seen, 
(i) Gauge invariance is satisfied, 

K';,q. = 0, (46) 

but spin conservation is not, 

9m^T- = - {l^S,o + ^q^5,^) i^AK^ ^ 0, (47) 

if r2 7^ 0; where i is a space component.— 

(ii) Depending on the relative magnitude of t^ and 
|Yct|, there are two regimes similarly to the charge chan- 
nel. More interestingly, however, for t^^ ^ |ycr|, the 
magnitudes of ps and jg can be independent, governed, 
respectively, by asymmetry in density of states and by 
asymmetry in conductivity; ps oc P^(7c and j'g oc (Tg. 

Finally, the spin density ps = {Js,o{q))uj and the spin- 
current density js = {js,i{q))u are given by 



-e2(A'" + A/v")div£;, 



= (7.E- 



{(iD^iyY) + 2T:f2{Du){(TDi 
Y^Y^ + 2T,f2{Yv) 



(48) 
-V(div£;). 
(49) 



2. Spin channel 

The response functions K^^^^{q,uj + iO) [Eq. (|A8p [ for 
spin density/currents are obtained as 



K: 



iO 



qMK' + AlC), 
q^LoK\ 



Kff = 



{aDv) ( % - ^ 



(38) 
(39) 
(40) 

^c.if«|,(41) 



3. Spin-resolved channel 

From Eqs. ([M]), ([S3, ([ISj) and ([H]), we obtain the 
"spin-resolved" density and current. 



where 



Pa = -e^K^dwE, 

j„ = cj^E + e'^DJ<„\I{AwE) 



D,Ys + 2tiT2{Dv) 

K„ ^ = Va 

Y^Y:^ + 2^T2{Yv) 



(50) 
(51) 



(52) 



From Eqs. ([50| and ([5T|) . we may derive 

> = (JaE - D^Wpa (53) 

where 



in the rotated frame, 



f,' = -eaV"o ( = JS). 



-e- — a)a°'Va. 



(59) 
(60) 



CTa = e DrrVrr 



(54) 



is the "spin-resolved" conductivity. Further discussion 
wiU be given in Sec. VI. 



The spin part of the impurity potential Vi^p is ex- 
pressed as Sficfa^c) = S^{t){a''a°'a), where Sf{t) = 

TZ"^{R'j,t)S'^ is the impurity spin in the rotated frame^i 
with 



IV. SPIN AND CHARGE TRANSPORT IN 
TIME-DEPENDENT SPIN TEXTURE 

In the previous section, we studied spin and charge 
transport in a ferromagnetic conductor in its uniformly 
magnetized state. In the second part of this paper, which 
consists of Sec. IV and Sec. V, we consider a more gen- 
eral case in which the magnetization varies in space and 
time. This magnetic texture and dynamics induce den- 
sity change and current even if A^™ is absent, which are 
calculated in this paper in the first order in both spatial 
gradient and time derivative. 



A. Transformation to local spin frame 

To treat the effects of space- and time-dependent mag- 
netization, we introduce a local spin frame where the spin 
quantization axis of s-electrons is taken to be the d-spin 
direction n{x) at each space-time point J^i^iii^ The orig- 
inal spinor c is then transformed to a spinor a in the 
new frame (rotated frame) as c = Ua, where L'^ is a 2 
X 2 unitary matrix satisfying c^(n • (t)c = a^a^a. It is 
convenient to take U = m ■ a with 



m 



sm - cos ( 



sm — sm ( 
2 



cos ■ 



(55) 



where 9 and (f) are ordinary spherical angles parametriz- 
ing n. From space/time derivatives, S^c = [/(9^ + iAp)a, 
there arises an SU(2) gauge field 



-iU^dM ^AT.a". 



ip — <.iv i^^i^ — ^1^,1^ . (56) 



This is an effective gauge field, which represents 
space/time variations of magnetization. The Lagrangian 
in the rotated frame is then given by L = Lei — H^-a, 



Ln 



Hn-f 



dr a' 



d h^ 

ih— + -— V^ + ep - Mmp + Ma' 
ot 2m 



dr ]ZAl 



2m 



dr AfAfa^ 



a'a, 



a, 
(57) 
(58) 



where j" = (/5", j") is a four current representing spin 
and spin-current densities ("paramagnetic" component) 



7e"^ = 2m"m'^ - ,5"^ 



(61) 



being a 3 x 3 orthogonal matrix representing the same 
rotation as U. Hereafter, the anisotropy axis of impurity 
spins is defined in reference to the rotated frame 



p^( Si ia,P = x,y) 



■^ ^ Si [a, 13 = z) 



B. Effective U(l) gauge field 



There is some arbitrariness in the choice of the rotated 
frame; one could take c = U'a' with U' = t/e"*'"'^/^, 
where x is an arbitrary function of x. This arbitrariness 
is a gauge degree of freedom in the sense that physical 
quantities should not depend on it. It is in fact expressed 
as the gauge transformation on a and A^, 



a' = e-^^^x/^a, 
A' = -i{U')^d,U' 



e^'^/^A^e' 



ia'x/'^ 



^"a^x/2, 



or, in componentwise. 



A': 



+ iA'y = 



{A 



. + ^K) 



K - a^x/2. 



(63) 
(64) 

(65) 
(66) 



Note that its z component A^ transforms like a gauge 
potential in ordinary electromagnetism, hence can be re- 
garded as a U(l) gauge field. In the following, when we 
refer to gauge transformation, it means Eqs. (|551l - ([55|) . In 
the next subsection, we study spin and charge transport 
driven by magnetization dynamics as a linear response 
to this effective gauge field ^^. 

Generally, one can do a gradient expansion in terms 
of A'^. The expansion parameter is qv^aTa and ujt„ 
(for A'l)^ where q~^ and lo are characteristic length 
and frequency, respectively, of the magnetic texture. In 
this work, we consider only the lowest nontrivial or- 
der in the expansion by assuming qv-p^T^ ^ 1 and 
wTg. ^ 1. This condition coincides with the condition, 
\Xa\ = \Dcq^ - iuj\T„ < 1, declared below Eq. ([T^. 
In typical experiments with Permalloy [v-pa ~ lO^m/s, 
T„ - lO-i^s)^, g-1 - 100 nm, w - 100 MHz^, we have 
Daq^T ~ 10~* and wTo- ^ 10~®, and the above conditions 
are satisfied quite well. 



C. Linear response to A™ and A^. 

Let us examine the density/current response to the 
two gauge fields, A™ and Af^. Spin density and currents 
considered here are the ones whose spin is projected on 
n (or z in the rotated frame), i.e., ps = p^ and js = j^ . 
The total current densities contain the gauge fields, 

.7, = {p,j + {epA'='^ + hp^A^)/m), (67) 

X,. = {ps/js + iepsA''''' + hpA')/m), (68) 

for charge and spin channels, where p = —ea^a and j = 

{—eh/2mi) a'^Wa. By generalizing Eqs. (fT8| and (fT9| . we 
may write 

{.]^.iq))^ = e^K^lA^ + eHK^lAl, (69) 

(JsAq))^ = e^A- Ar + efti^- AJ. (70) 

The response functions, K"l and A'^"^, are obtained from 
Eqs. (PO)) and (PT|) by replacing the electron operators in 
the original frame, c (c^), by those in the rotated frame, 
a (a^), and are already calculated as K"^^ and A'^'J, in 
Sec. in. Thus the response to ^™ in Eqs. ^ and dTD]) 
exactly follows the results there. 

Let us then focus on the response to Af^, in particu- 
lar, on ATJ;^. (ATj^'J^ will be presented in Appendix D.) 
From the definition (linear-response formula), one can 
show that the Onsager's reciprocity relations hold. 



K';^{q,tUJ^) = K\ 



lyfi 



{-q,-iujx), 



Kf:^{q, LO + iO) = ATJ-q, -lo - iO). 



From this, we see that 



g^if^t = KlU^^ = 0, 



(71) 



(72) 



(73) 



namely, the charge conservation is satisfied also in the 
response to A^^. On the other hand, if r2 ^ 0, spin 

is not conserved, q^K^^ 7^ as seen before. This fact, 
combined with Eq. ([7^ . implies that A'Jj^ is not gauge 
invariant, 



if- g, = 9,X- ^ 0, 



(74) 



if r2 7^ 0. The gauge non- invariant terms in Eq. ([55)1 
may be extracted as^ 

j;(q,w) = eh^K^{q^5^o + g.wJ^JA^^Q. (75) 

To summarize, the calculation based on the gauge field 
A^ fails to respect gauge invariance in the presence of 
spin-flip scattering. Stated more explicitly, the density 
and current calculated as a linear response to A^ are not 
gauge invariant 1^ 



V. CAREFUL TREATMENT OF SPIN 
RELAXATION EFFECTS 

A. Restoration of gauge invariance 

The lack of gauge invariance encountered in Sec. IV- 
C is due to an oversight of some contributions. We re- 
call that the quenched magnetic impurities in the origi- 
nal frame become time-dependent in the rotated frame, 
Sjit) = TZ"^ {R'j,t)S^. Therefore, we should treat the 
spin part of the impurity potential 



H,=us^f drSj{t)5{r - R'^) ■ (atcra)^ 



(76) 



as a time-dependent perturbation. The same situation 
was met in the calculation of Gilbert dampingJ^ 

Since the first-order (linear) response vanishes, 

S°'{t) = 0, let us consider the second-order (nonlinear) 
response. 



Aj^(g,'^) 



tgLtg 



duj' 



-ensu: I —xfiq;i^,i^')[SH^ - ^')SH^')], 



(77) 

where Sp{uj) is the Fourier component of ^ ^^"(i) S{r — 

R'j), and Xn^ is the nonlinear response function.— To 
calculate it, it is simpler to use the path-ordered Green's 
functioui^ The contribution represented in Fig. [5] are 
given by 

xtr[(tv + A^)Gfe+(e+)a"Gfc.(£ + ^')^''Gfc_(e)]< 

(78) 

where e+ = e + uj. The Green's function Gk{e) now 
stands for a path-ordered one, whose lesser component is 
given by 

G<ie)^fie){Gtie)^G^ie)), (79) 

with /(e) being the Fermi distribution function. In 
Eq. ([75)l . we adopt a matrix notaion, {G)ccr' ~ GcSca', 
{K)a,a' = AJAa' with A;^ given by Eq. |M|, and 'tr' 
means trace in spin space. 

We expand Xu^{Q't^,^') with respect to u and to' as 



xf (g;c.,c.') = Af - lioBf - iLo'Cf + ■ 



(80) 



where A"f , B"^^ and C"^ are the expansion coefficients. 
Substituting Eq. dHHl) into Eq. ([771), we have 



Bfdt{S"SI3) + Cf S°'dtSP 



(81) 



xf = ^^ 




e + u) . a 



e + uj' 




£ + U)' 



FIG. 2: Diagrammatic expression of Xm ■ The wavy line 
represents scattering from impurity spins, which are time- 
dependent in the rotated frame. The shaded part represents 
the vertex function A^. 



where S = S{t) is time dependent. (We have dropped 
a term containing A"*^, which does not reflect the time 
dependence of S{t).) From 



S'^dtSP = {Sl5°[^ + SI 5°"'5^''){ndtny'', (82) 
where (5" = 5°"^ — 6"^S^^, and the relationii 

(TZd^n)''^ ^2e''^^Al, (83) 

we see that Eq. (|8T|) describes a response to Aq. The 
coefficients are calculated as^ (see Appendix B) 



= ,,.,, i^lH,i±M^:^^^e«^ (84) 
Y^Yi + 2t:T2{Yi') 

where e"" = e"^^, and we have dropped unimportant 
terms proportional to (5" or (5"^ (5'''^. We thus have 



Aj^iq^io) = ehAK';,Al, 



(85) 



with^ 



Ak^l = -AK^qHf,o + q^i^S^^}S,o■ (86) 

This new contribution cancels the gauge- dependent 
terms, Eq. (|75p . and restores the gauge invariance, 



(A-t + AA-t)Q. = 0. 



(87) 



Note that it does not affect the charge conservation since 
q^AKf^l, ~ 0, nor the spin non-conservation {q^Kfl^^, ^ 0) 

since it does not contribute to ATi^'^,. 

The gauge-invariant result for the charge density 



{q,u}) and current density j''™^(^)(q', w) induced 



psmf(i; 

by magnetization dynamics is summarized as 



•smf(l) ^ _gO 



+ ,.(^^!^^i>±Mi(M(^Mv(div^O). 

Y^Yl + 2TTr2{Yi^) 

(89) 










+ 1^, 



+ aj' 



CD. 



+ V, f I) + ^1' 






a 

-Lo a" 




e + u' 



a 'Vi a' 



ji 




FIG. 3: Diagrammatic expression of xZi ■ The gray circle 
represents the interaction with AJ^. 



The first term on the right-hand side of Eq. (|89[) has 
the form of Eq. ([3]), and implies the existence of spin- 
dependent motive force described by the effective 'elec- 
tric' field E^. The second term of Eq. ([5^)) represents a 
diffusion current arising from charge imbalance induced 
by £J°, as made clear in Sec. VI. This term implies the ex- 
istence of nonlocal spin-transfer torque as the reciprocal 
effect, whose study will be left to the future. 



B. Dissipative correction 

It is important to note that there is one more con- 
tribution within the same order in gradient expansion. 
It is essentially given by Eq. ([77| . but with one more 



factor of A". The response function, denoted by xZi ^ 
is obtained from Eq. ([75)) by further extracting A" via 
Eq. ([55]) . These are expressed as [Fig. [3] 



3r'^'\qM 



-ehriM, 



2 Y^ 



s"'s /_^ 






X [S-iu;~u')Sf^icv')] Al 



where 



x;'^-(g;c.,c.') = E 



k,k' 



de 
27ri 



tr 



q-q' <?',«' 



(«M + K) 



(90) 



X {v+Gk+{e+)a^Gk+{e+)a''Gk'{e + io')<j"Gk-{e) 
+vrGk+ (£+)a"Gfc. (e + io')a^Gk- {e)a'^Gk- (e)} 

+ - 5^,a-'Gk+{e+)cj^Gk,{£ + Lo')aPGk-{e)Y . (91) 
m J 

with vf = {ki±qi/2)/m. We have put q' = in Eq. ([9T|) . 
but retained q and uj. Note that the terms with 7 = 2 
cancel out, and A^ does not contribute. In the same way 

as Sec. V-A, we expand xZi with respect to uj and lu' 
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as xfi^ ^ Alp - iujB'^f^ - i^'C^f'' + ■■■ and focus on 

the coefficients i?" and C" . Deferring the details to 
Appendix C, we cite the result 

1 



a/37 ^ j^iCt/S'y 

(5-e^7„^/3.^a7)i^^^(i^^)RA^ (92) 






where j/+ = i/^ + vi, and Lf„'s are given by Eqs. (|A15P 
and (jA17p . Note the order of the subscripts, ifi. We thus 
have 

jrf(2)(q,c.) = /3^^a(Lr,f^(A^.Ao^),.„ (93) 



where Aj^ = A^ — z (z- A^J, and 

/3=^"sHs'(^ + ^)(^t + n) 

is a measure of spin relaxation. With the relation 

1 



-^,i ' -^n 



■ h- diTi^ 



(94) 



(95) 



which is gauge-invariant under (|65l) . we finally obtain 



psmi 


(2) 


= 


-e2i^sdiv^f^ 




(96 


jsmf 


(2) 


= 












+ 


2 {aD^iyY) + 2TTt2{Du){(7Dv) 


V(diA 


pdisN 




Y-^Y^ + 2tiT2{Yv) 


f^a J' 



(97) 

where E'^^^ is given by Eq. ([2]). Since E^^^ contains j5 as 
a prefactor, Eqs. ((96|) and ([97)) come from spin-relaxation 
processes. This /3 is exactly the same as the coefficient 
of the /3-term of current-induced torquCj^iii consistent 
with the fact that these are reciprocal to each other »^i^ 



VI. RESULTS AND DISCUSSION 

The results obtained in this paper are summarized as 

{DvYF) + 2'kV2{v){DvF) 



j = acE + c7sEs 

{D'^vYVF) + 2Tif2{Dv){DvVF) 
^ Y^Y^ + 2^f2{Yv) ' 

{(jDvYF) + 27rf 2 {<7v) {DvF) 
^' ~ Y^Y^ + 2^f2{Yv) ' 

js = CTs-E + (TcEs 

{aD^vY^F) + 2Tif2{(TDv){DvyF) 



+ 



Y^Y^ + 2^T2{Yu) 



(98) 

(99) 
(100) 

(101) 



where F^ = e^div.E, Fg = e^div.Es, and F^ ~ Fc + aFg. 
The notations are as before; for example, {DvVF) = 
D^^v-^VF-^+D^v^VFi. From these relations [or Eqs. ([T04[) 
and p05p below[, we identify the spin motive field to be 

E,,, = — {-n-{ri X d,n) + /?(n-9,n)} . (102) 
Ze 

The spin-resolved density and current are given by 

p, = -e2div(A'^£; + A^;£;,), (103) 

> = cj^E^-D^Vpa, (104) 

E„ = E + aEa, (105) 

where E^ is the total field felt by spin-a electrons. The 
coefficient K^, is given by Eq. ([52[) . and K'^ by 

K = ^^.Z£±^^£a(^M.^. (106) 

Y^Y^ + 2^T2{Yv) 

There are two characteristic regimes depending on the 
relative magnitude of r^j^ = 2'KV2^v)lh and \Xa\- For 
Tg7^ < IF^I, Eq. pl?5)) becomes 



P'T 



i?.?^ 



■ div.E;„ 



(107) 



meaning that the spin- cr electrons respond only to E„^ 
not to Eg^ and the two spin components (^ and 4^) be- 
have independently. In particular, the response to a spin 
motive field E^ (set E ~ ^ for simplicity) is opposite in 
sign between ^ and ^ electrons. In the opposite limit, 
tJ^ > jy^l, Eq. pl?g| becomes 



v^l{v) 



D,m^ 



■d:w{a^E^ + a^E^), (108) 



where D^s = {Dv) / {v). In this case, the density of spin- 
a electrons is affected not only by E^ but also by E^. 
This is due to the strong spin mixing; as an elementary 
process, p^ is induced solely by E„, not Eg^ but subse- 
quent spin-flip processes tends to equilibrate p^ and p^. 
Note that f electrons and \. electrons respond to Eg with 
the same sign. (The common sign is determined by that 
of (Tf- -(7;.) 

The above features oppose the picture of two indepen- 
dent currents, but they are actually described within the 
conventional two-current model.—"— This is best demon- 
strated by the relation 



di 



Pa + divjcr 



where 



'si, a 



Pa 

Tai,a 



2T:T2l^&/h 



_Ps_ 

Tai,S 



(109) 



(110) 



is the spin-flip rate for spin-cr electrons. The right- 
hand side of Eq. p09p represents a coupling between f 
and \, electrons. In deriving Eq. (|109p . we have used 



Eqs. (UnSl), (PM)1 . ([TU5)) and ^, and the relations, 
(ctA'/i/) = {aDY)/Z and {(tK'/v) = {DY)/Z. Note that 
Per, being given by Eq. (|103p . represents a deviation from 
the equihbrium value. One may define the deviation of 
chemical potential, 5^,^^ from equilibrium by 

Pa = -ev„5pL„. (Ill) 

Then Eq. (|109p can be put in a familiar form^i^ 



d , J. . O'er 5p,a - S^^ 



(112) 



where ^o- = y^TJyr^i^ is the spin diffusion length for spin- 
a electrons. 

The present work is therefore within the two-current 
picture. This fact was implicitly used in identifying the 
spin motive force on the basis of Eq. ([3]). 



VII. SUMMARY 

In this paper, we have studied spin and charge trans- 
port in a conducting ferromagnet driven by two kinds of 
gauge fields, A™ and Af^, which act in charge channel 
and spin channel, respectively. In particular, we have 
given a microscopic calculation of spin motive force by 
taking spin-relaxation effects into account. 

In the first part, we calculated density and current in 
both spin and charge channels in response to the ordi- 
nary electromagnetic field A™ in a uniformly magne- 
tized state. We observed a crossover from two diffusion 
modes to a single mode as the spin-flip rate is increased 
(for a fixed frequency/wavenumber of the disturbance), 
or as the frequency/wavenumber is decreased (for a fixed 
spin- flip rate). However, if expressed in terms of spin- 
resolved density and current, the so-called two-current 
model is shown to hold irrespective of the strength of 
spin-flip scattering. 

In the second part, we have developed a microscopic 
theory of spin motive force in the framework of gauge- 
field method. We readily encountered the problem of 
gauge non-invariance; the current calculated as a linear 
response to A^ depends on the gauge (choice of local 
spin frame). This fact is intimately related to the non- 
conservation of spin (due to spin-flip scattering) by On- 
sager reciprocity, hence is robust. This theoretical puzzle 
was resolved by noting the fact that the spin-dependent 
scattering terms (quenched impurity spins) are time- 
dependent in the rotated frame. By calculating the 
second-order (nonlinear) response to this time-dependent 
perturbation, we could recover a gauge-invariant result 
while keeping the spin non-conservation. The dissipative 
correction to the ordinary spin motive force, which is the 
inverse to the spin-torque /3-term, is also obtained. 

Note added: After submitting the manuscript, we be- 
came aware of a closely related work by Kim et al^ 
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Appendix A: Calculation of response functions K'^^^ 

and Af,";, 



In this Appendix, we evaluate the electromagntic re- 
sponse functions in the ladder approximation shown in 
Fig. m^a). From Eqs. ((20|) and ((2T|l . they are written as 

Kll{q,iwx) = -TY,L%{q-Aen + i^x.ien), (Al) 

n,<7 



with 



-^5^1/ — nj^i/ + nj^o^S^' 



(A3) 



^% =^'"tJ.^>^Gk+,aii£n+i^\)Gk_,a{i£n),{J^'i) 



where e„ = {2n+l)TTT {n: integer) is a fermionic Matsub- 
ara frequency. The vertex function A'^ satisfies [Fig.[T]^b)] 



A^-A^ + Tin^A^ + rsn^A^, 



(A5) 



where 11^ = II^Jq, and AJ^ == Tili;^^ + r2n^^ is the lowest- 
order contribution. The equation (|A5[) is solved as 



A'^ = 



A^-n^(riA^-r2Ag) 
1 - fi(nt + H;) + (f2 - f2)ntn^ ' 



(A6) 



Performing the analytic continuation, iuj\ -^ uj + iO and 
retaining terms up to the first order in uj, we obtain 

X-(q,w + zO) = ,.+<5^o^,o + ^^(i-jRA^ (A7) 

(7 

lC;,{q,Lo + iO) = ^_5^o5„„ + ^^a(L^J^A,(A8) 



where iy± = i/^ ± v^. The function (L'^^)^^ is obtained 
via the analytic continuation, z(e„ + u!\) — >■ s + uj + iO 
and ie„ ^- e — iO, as indicated by the superscript "RA". 
We assume j^ ^ £Fa, and discard (LJ^^)^^ and [L'l^)^^ 
as in usual calculations of transport coefficients. The k- 
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integrals are evaluated up to OdXcrl) or ©(iXo-l^) as and thus 



k 






fc 

~ 27ri?,:/^%. (All) 

where D„ = Up^ro-/3, wpo- = hk^cr/in, and Xo- = YctTct 

with Yo- = D^q^ — iuj. Using these formulas, we obtain t^ , n _p ; 



rt>4, + 27rr2(ri^) 



-tut- 



RA_ r. + 2^f,M 1 ,,^^, g^ yrn + 27rr2(y^> 

"'' yt^4. + 27rf2(rzy) T^' 



(A17) 



,,^_^ D,lW:2^rMM.l, (A13) 

yt^4, + 27rr2(y!^) t^ 



Appendix B: Calculation of C"*^ 

The nonlinear response function Xa^ in Eq. (|78p is written as 






— (i^^(q;£ + c., £)/.,(£ + c.'))<, (Bl) 



where Lg is given by Eq. (jA3|) . and /o-(£) = X^fc Gfco-Ce)- Following the Langreth's niethodr^i^ the lesser component 

of Lq (q; e + Li!, e)Ig-{e + uj') = LI is calculated as 



{LI)< = /(£)(i^^ - L««)/R + /(e + uj')L^^{I^ - /«) + /(£ + w)(L^^ - L«^)/^. (B2) 

Note that the ordering of Green's functions in LI is G{e + w)G'(£ + w')G(£) [see Eq. ([75)) ]. The superscripts RA, A 
etc. specify the analytic branch; for example, L^^{e + uj,e) = L{e + w + iO, £ — iO), /^(£) = I{e — iO), etc. Thus the 
coefficients in the expansion x'ff ~ ^'ff ~ iujB"^ — ito'C^^ + • • • are obtained as 

(T.CJ' 

Cf - -E[(<5f + za£«'')<5,,^ + <5"^<5^^J.,,](L^^(q;c.,0))«'Alm/«,(0). (B4) 



tr.cr 



We have retained only the lowest-order term in 7o-. Substituting Eqs. (|Al4p and (JA16P together with /^(O) = —i-Kv^ 
(whose real part is dropped consistently with the selfenergy) into Eq. (jB4[) , we obtain Eq. 
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Appendix C: Calculation of C" 

Consider the nonlinear response function x^ii given by Eq. (|9ip . As in Appendix B, we take a lesser component, 
extract the w'-linear term, and retain terms containing both G^ and G^ to obtain B" = — (l/2)C"f''' and 






c. -^J2tr[{v, + A«^)«,{G«+a^G«+a"i>a^Gt + G^+cT"t>a^Gta^Gt}] 

k 

~^S,^Y,tr[<J''G^^a"OaPGt]■ (CI) 

Here {A^^)„„, = (A;^)^^(5^^, is given by Eqs. (Pl2| - ((IT3|l . and i> = Efe'(Gfc. - Gg,)/27ri is a matrix of density of 
states, {i>)cra' ~ I'aSaa'- In Eq- (|Cip . all G's are evaluated at e = except for those in A^ in which q,Lu are retained. 
Equation (JC1|) is written as 



where 



M-(q,..) = Q;,iq) + iK)'''^Qo^iQ), (C3) 

Af;,(g,c.) = (3^,(g) + (A^)«'AQ[;,(q), (C4) 



In the lowest order in 'y^, we see that 

M;,{q,u) = M;,{q,u) = -^m^)^^, (C6) 

where (Lf^)^"^ is given by Eqs. (|Il5)) and ((JTT)) . Noting that A/;^^ + ii/*^ = - E^ ^{^1^.)^'^ /"^^'i is independent of 
(7, we obtain the leading term as 
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Appendix D: Spin current induced by spin motive Hence, we have 

force 



The response function K't^ in Eq. ([70]) is evaluated as (j^a- nRA ^ 2nv °"'^°' ^ '^^^^ ^'^^' fD6) 



is"u. = am,+n- A- . (D2) 



^<x xRA . . (jD,Ys + 2TTr2{<jDl^} 



{L^,o^) = -2^iq^u, ^^ .r^ '-, (D8) 

The spin-current vertex function A^^^ , which satisfies ^t^i + 27rr2 (Yv) 

K. = K. + f in.A^,, - f 2n,Af,,, (D3) (x^,^.)«a ^ sttT^.i.. L<5,, - q^q^^£^^±Mlj^\ . 

[ y^y; + 27rr2(ri^> J 

with A^,, = fT(fin[J, - fsHgJ, is given by (D9) 

(Xa ^RA _ g^ff + 27rf 2 (cTi^) 1 

l^'-s.oJ — vva-9r/v\T' l^*/ Note that K'f^'s thus obtained do not satisfy spin cou- 



rt >;+ 27rr2 (Fi^) Ta 



servation nor gauge invariance, QfiK^j^ = K^^q^ ^ 0, if 



Y^Y\^ + 27rr2(y!^) To- 

I 

As in Sec. V, time-dependent magnetic impurities, Eq. (|76p . in the rotated frame also induce a spin current 



Ajs,^(g,^) = -en^w,' / — x"2(g;^,^')[^"(^-^')^^(^')], 



where 



/OO 1 / 

^ x"^7(g; c., c^') [^"('^ - c^')s^{co')] ^.A^,,^, (Dio) 

-OO '^^^ 

^ tr[Ka^ + As.^)Gfe^ (e + a;)a"Gfc. (e + u:')a^Gk_ {e)]< , (Dll) 

/OO I 

-^ tr [(z;^a^ + A,.^)z;+ G^^ (e + uj)<j-'Gk^ {e + u)a''Gk' {e + iu')a^Gk_ (e)] < 

K,K '°° ^* 

/OO I 

— tr [{v,,a' + A3,^)z;- Gk+ (e + oj)a"Gfc. (e + io')a^Gk_ {e)a'<Gk_ (e)] < 

-I /'OO J 

+ -5^,P^Y. tr[Gfc^(£ + c^)a"Gfc-(£ + c^'y'^Gfc_ (£)]<, (D12) 



As,/xz 



with Wj = (fci ± qi/2)/m. We have put g' = in The first term in Eq. (JDISP corrects (the first two of) 

Eq. (|D12p . By taking the lesser component and extract- the following response functions, 
ing the w- and oj'-linear terms, we have 

(D13) 

with kll + ^K^Q = q'^Ki, (D15) 

{Y)5,o-^q.{DY)S,. , ^^o + A^S = MiDi^) l'^,}, (D16) 

AK^^ = -4.r2.,n y^^^^2.f2(r.) '''°- i^;;? = ..-A, (D17) 

(D14) i^^l = iuj{{Diy)S,,-q,qjK2}, (D18) 
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where 



K, = (^-n^+2.r.M(^i^^)^ (019) 



Ko 



Y^Y^ + 27rt2{Yiy) ' 



(D20) 



and restores the gauge invariance. This leads to a spin- 
current density, 



-fW(q,..) = ^Y.^{Ll,^r'^El (D21) 



The second term in Eq. (JD13I) gives 



-f(2)(g^^) ^ |^^^(i^^^^)RA^dis^ P22) 



Therefore, the total spin-current density induced by the 
total spin motive field Eg = E^ + E^^^ is given by 



•sinf 



(D23) 

jT' - a,Es + e'K2VidWE,). (D24) 
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